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TECHNICAL MEMORANDUM 


MEASUREMENTS AND CALCULATIONS OF THE COULOMB CROSS SECTION FOR THE 
PRODUCTION OF DIRECT ELECTRON PAIRS BY ENERGETIC HEAVY NUCLEI IN NUCLEAR 
TRACK EMULSION - CENTER DIRECTOR'S DISCRETIONARY FUND FINAL REPORT 


I. INTRODUCTION 

The principal goal of this study has been to develop a practical method 
which employs only the electromagnetic interaction to measure the energy of 
ultrarelativistic heavy ions in the energy region above 10 13 eV/amu (10 
TeV/amu). The only apparent electromagnetic phenomenon that may be used is 
the materialization of electron-positron pairs from the virtual photon field 
of the relativistic ion in the field of a stationary nucleus. This is a 
phenomenon which has been studied experimentally for more than 35 years 
without resolving differences between measurements and calculations, and thus 
is still an interesting aspect of quantum electrodynamics ( QED ) • Until the 
last several years, work on this phenomenon has been limited to primary 
particles from cosmic ray air showers (electrons and mu-mesons), a few heavy 
primary cosmic ray nuclei, and protons and TT-mesons from accelerators. The 
yield of direct Coulomb pairs from the singly charged primary particles has 
been less than one per meter in the typical production and detection medium 
(nuclear track emulsions), and thus a definitive experimental study has been 
hampered. The recent availability of heavy ions with high energy at the 
European Center for Nuclear Research (CERN) has opened the opportunity to 
resolve the previous experimental uncertainty and confirm the yield of direct 
Coulomb pairs. 

Present techniques for measuring particle energy which employ 
electromagnetic phenomena such as Cerenkov counters and the relativistic rise 
of ionization in gas-filled counters (ion chambers and scintillators) saturate 
below 1 TeV/amu for practical detectors. Transition radiation detectors have 
been designed to measure cosmic ray energies to about 10 TeV/amu. Only 
"ionization calorimeters" that depend on the strong (nuclear) interaction have 
been accepted as practical energy measuring devices above 10 TeV/amu. 
Calorimeters require enough depth that the primary particle interacts 
(nuclear) in the device and produces a cascade of hadronic (p, tt, K, etc.) and 
"electromagnetic" particles (e, y), the number of which are sampled to 
estimate the primary energy. At 10 14 eV/amu this requires at least 5 cm of 
lead for thin lead-x-ray film-emulsion calorimeters, >30 cm of lead for "total 
absorption" calorimeters of lead-scintillator, or >15 m water equivalent depth 
in the atmosphere for air shower arrays. In addition to the mass of material 
required for ionization calorimeters, there are serious concerns about 
possible changes in the behavior of nuclear interactions beyond 10 14 eV, 
particularly for heavy nuclei. Part of these concerns has been generated by 
conflicting results of air shower experiments which, by various analysis 
techniques, have indicated that cosmic rays are either mostly protons or 
mostly iron near 10^ eV. Other concerns about nuclear interactions at high 
energy stem from the theory of strong interactions, quantum-chromodynamics 

(QCD). QCD predicts for heavy nucleus interactions a phase change of nuclear 

matter from bound quarks in neutrons and protons to a plasma of the 

constituent quarks and gluons. This would probably result in a quite 

different particle cascade in calorimeters. 



A technique for energy measurement that does not depend upon nuclear 
interactions would be most desirable in the energy regime above 1 0 1 ** eV/amu. 
The potential technique we have investigated is the measurement of the linear 
frequency of occurrence of direct Coulomb pairs along the track of heavy ions 
in emulsion [1]. Calculations described below indicate that cosmic ray iron 
nuclei at 10'** eV/amu should produce ~30 direct pairs per centimeter of 
pathlength in nuclear track emulsion. A pathlength of 3 cm of emulsion shouLd 
allow an energy measurement to about 50 percent accuracy, which is competitive 
with the calorimeter technique. However, most previous measurements of the 
electromagnetic pair yield reported in the literature were not consistent with 
calculations and the past calculations have involved questionable assumptions 
and approximations. The calculations generally converge at the highest 
energies, but give different results below ~1 0 eV/amu. Consequently, we 
have approached the problem both with experiments and by examination of past 
theoretical approaches and improved calculations. The experiments have 
involved exposure of emulsions to heavy ions at particle accelerators with 
various energies. Iron ions at 2 GeV/amu (at Bevalac, Berkeley), oxygen at 60 
and 200 GeV/amu (at the CERN Super Proton Synchrotron, SPS, Geneva), and 
sulfur at 200 GeV/amu (CERN) have been used to produce >100 meters of ion 
tracks for each exposure. Much of this has been scanned with microscopes to 
produce results as described below. 

In the next section the results of previous experimental studies in the 
literature are summarized. These results are compared with our experimental 
results which have recently increased in pair frequency as a result of 
"learning experiences." These experiences have involved both scanning 
technique and attention to the results of theory concerning the energy and 
angular distribution of the electrons in the pairs. 

The scanning procedure described in section IV appears adequate to 
efficiently find pairs down to a total pair energy of about 3 MeV. This is 
sufficient to establish a technique, provided a long pathlength in emulsion 
(>1 cm) is available and the dip angle of the track is not too steep. 

The principal problems that remain in establishing the pair method for 
measuring energy are the uncertainty in the calculation of the yield for the 
low electron pair energies (<10 MeV) and concern whether the latest 
measurements in this study are contaminated by any significant background 
effect. The calculations in this study are summarized below and reported in 
detail in published references. Further work is planned to attempt solutions 
from exact QED to resolve present ambiguities for the low electron pair 
energies. 


II. SUMMARY OF EXPERIMENTAL RESULTS 

The direct production of electron pairs by energetic ions in nuclear 
track emulsion was first studied experimentally in the 1950's [2-5]. 
Relativistic electrons in showers of particles caused by cosmic rays served as 
the primary particles, whose energy ranged from lO^-IO 5 MeV (see Figure 1). 

The difficulty in analyzing these "Trident" events [e_ -> e_ + (e_, e + )] was 
due to a high frequency of bremsstrahlung photon production by the primary 
electron and the conversion of those photons to pairs [e_ -> e_ + ( y -> (e_, 
e + ))], which required a large correction. Thus heavier projectiles are more 
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appropriate primary ions since the production of bremsstrahlung will be 
reduced according to the ratio [m e /M] 2 . Energetic muons, pions, and protons 
were studied [6-11] after particle accelerators, which produced high energy 
monoenergetic beams, became available in the late 1960's and early 1970's. A 
beam of 200 GeV protons at the Fermi National Accelerator Laboratory was 
utilized to produce direct electron pairs in an emulsion [10-11]. The results 
are included in Figure 1. A small sample of heavy cosmic ray tracks, taken 
from emulsions exposed on "Intercosmos^" [12], was examined for direct pairs 
and the yield is included in Figure 1. This cosmic ray data of Grigorov et 
al. [12], although limited by poor statistics, inexplicably reported a higher 
value than the proton accelerator measurements even though the reported pair 
energies are very high. Fortney et al. [9] used a Ne-H 2 bubble chamber in a 
25-kilogauss magnetic field to measure the momentum and yield of pairs above 
10 MeV/c. The results matched their calculations of the yield and pair 
momentum spectrum. If adjusted for the high pair momentum cutoff, and for the 
target and the projectile atomic numbers, this result would fall close 
to (A) in Figure 1. However the appropriate pair momentum cutoff adjustment 
is presently uncertain as discussed in section VI. 


Heavy ion beams have become available in recent years first at 2 GeV/ amu 
(Berkeley), then at ~10 GeV/ amu (Brookhaven ) , and more recently at 60 and 200 
GeV/amu (CERN-Geneva ) . We exposed nuclear track emulsions to all these beams 
and have scanned significant track lengths of most exposures to determine the 
linear frequency of pairs. Results of various individual measurements in this 
study are included in Figure 1. Present direct electron pair yields have 
increased commensurate with the ability to recognize the lower energy pairs. 
The initial problem in detecting these low-energy electron pairs (e < 10 MeV) 
was due to the fact that the division of the energy between the electron and 
the positron can be so disparate that it was difficult to recognize the lower 
energy electron (positron) as a member of a pair and to trace the lower energy 
electron (positron) to a common interaction vertex. Although the lower energy 
pairs are generally accompanied by larger angles between the primary particle 
and the members of the pair they are not strictly correlated to pair energy, a 
consequence of the five-body nature of the physical phenomenon. These 
aspects, covered in detail below, make the lowest energy pairs difficult to 
recognize. This is especially true of electrons (e_ < 1 MeV) that can either 
be mistaken for delta rays or are missed altogether due to the background. 

The data from proton tracks cited above are consistent with this assessment 
since they detected only the highest energy pairs resulting in a low yield. 

The recognition of these problems resulted from new calculations of pair 
energy and angle distribution and also a new scanning technique and electron 
energy measurements described below. 


III. SUMMARY OF CALCULATIONS 

To better understand the physical properties of the direct electron 
pairs, calculations have been performed that predict the total yield, the 
energy transfer during the collision, the energy sharing, the net transverse 
momentum, and the angular distribution of the pairs. The standard theory is 
defined to be the lowest order quantum electrodynamic (QED) evaluation of the 
total cross section for pair production by heavy ions in emulsion. Calcula- 
tions [13-20] dating back to the 1930's include both the classical and QED 
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treatments of this problem where some approaches have taken into account the 
atomic screening of the absorber. Computational efforts reported here have 
focused on the Weizsacker-Williams (WW) method of virtual quanta as described 
in [21], including corrections for pair production by low-energy photons 
(Racah formula [22]) and higher order corrections [23,24] for relativistic 
pairs (i.e., e + >> 1 and e_ >> 1). Also included are higher order terms 
involving the effective atomic number of the target. The region of validity 
for the WW method is determined by two conditions: (1) the Lorentz factor for 
the heavy ion (y) is much greater than one, and (2) the energy transfer 
(e/m e c 2 ) is much less than (y). One of the shortcomings of the WW method is 
that it depends on an undetermined parameter which corresponds, in this case, 
to the minimum impact parameter. Fortunately another method that evaluates 
the direct pair production cross section (Kelner and Kotov (KK) [18]) has no 
undetermined parameters but is valid only in the relativistic energy transfer 
region. If one matches the results of the KK and the WW calculations in an 
energy transfer region where both should be valid, then the minimum impact 
parameter is determined (see P. Eby [25] for more details). The total 
estimated yield for y = 200 is shown in Figure 1 (+ symbol) and should be 
compared with the most recent measured yield (sulfur - solid triangle) that 
reflects the current scanning criteria. The corresponding total energy 
transfer curve (dN/dE) is compared to the same set of electron pair data taken 
from the CERN sulfur exposure (Figure 2). It is quite evident that there is a 
measured overabundance of the low-energy pairs compared to the theoretical 
estimate. This discrepancy is being examined from two sides. Do the 
computations accurately predict the frequency of the low-energy pairs and is 
there any systematic bias in identifying the low-energy pairs in the 
emulsion? These questions are discussed in section VI. 


IV. DESCRIPTION OF ACCELERATOR EXPOSURES 

The EMU04 emulsion experiment design [26] is illustrated in Figure 3. 

The basic configuration is a stack of 20 thick emulsion plates (800 y acrylic 
substrate with 500 y of Fuji 7B emulsion coated on both sides) oriented 
parallel to the beam and housed in a light-tight black plastic box with a thin 
window facing the beam. The emulsion pouring was accomplished at the Fermi 
National Accelerator Laboratory in Batavia, Illinois. The total exposure of 
the emulsion plates to the heavy ion beam at the CERN SPS was controlled by a 
programmable deflector magnet. The magnet was positioned approximately 2,000 
m upstream of the target area where a preselected fraction of the heavy ions 
are deflected from the main beam. The normal beam intensity had to be reduced 
so that the ion events could be analyzed in the emulsion without interference 
from neighboring events. This was achieved by defocusing the beam so that a 
beam intensity of 2,000 to 3,000 ions cm" 2 was obtained. Thus the central 
intensity, lateral spread, and ion purity were carefully monitored. Figure 4 
shows the beam access area where each of the EMU04 emulsion stacks were 
exposed to two separate beam spills displaced 1 cm apart. The TV monitor 
displayed the beam profile which assisted the operator in timing the insertion 
of the emulsion targets. Just after the exposure and before the plates were 
developed, a reference grid was photographed on one side of each emulsion 
plate. This grid contained a pattern of numbers that allowed unique 
identification of primary tracks. 
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Over the past several years the EMU04 emulsion plates have been exposed 
to the 60 and 200 GeV/amu oxygen beams and the 200 GeV/amu sulfur beam at the 
CERN SPS • Preceding the CERN heavy ion exposures, emulsion stacks of a 
similar design were exposed to a 1.8 GeV/amu Fe beam at the Lawrence Berkeley 
Laboratory (LBL). Figure 5 is a microphotograph mosaic of a typical direct 
electron pair in the emulsion oriented close to the focal plane of the 
microscope. The heavy primary ion is S 32 . 


V. MICROSCOPE SCANNING: TECHNIQUES AND RESULTS 


A total of 50 m of heavy ion tracks have been scanned with a 
magnification of 1000. Figure 6 shows the measured pair yield as a function 
of distance from the edge of the emulsion plate. The yield is being gradually 
attenuated as the track length increases. This is caused by (1) a nuclear 
interaction prematurely terminating the direct pair search, (2) the background 
level building up with track length in the emulsion making it more difficult 
to observe the pairs, and (3) the tracks not being perfectly parallel with the 
substrate (i.e., the dip angle is not zero); therefore, the lateral spread 
results in a distribution of track lengths. Hence the first scanning rule is 
that the pathlength/track is limited to 3 cm from the edge of the emulsion 
plate. The initial set of measurements of pair frequency in this study is 
shown in Table 1. These results, normalized to 3 cm of track length and to Fe 
ions (by Z 2 ), are plotted in Figure 1 (A,0,0,(j)). They are above previous 
measurements with proton beams, but below previous and our recent 
calculations. Examination of calculations of production angles and energy 
sharing between electrons indicated that low-energy electron pairs were being 
missed. An exercise with different scan rates (shown in Table 1) indicated a 
scan rate of ~1 cm/hr should be the upper limit with experienced scanners. 
Because of the low yields of these measurements a new scanning procedure was 
developed by D. T. King [27]. Additionally, the measurement of the energy of 
each electron was adopted, as well as their angles with respect to the 
primary. The scan procedure and electron energy measurement procedure are 
described below. The result of scanning with this procedure is given in 
Figure 1 (*). 

The latest procedure is as follows. Scan at 1000 magnification at a rate 
of <1 cm/hr, focussing at each field of view over the depth >20 y above and 
below the primary particles track. Note every knock-on electron of energy >2 
MeV (the range is at least one field of view and in the same general direction 
although it may be scattered). The energy of the electron is measured by 
recording three deflection readings for the multiple-scattering angle over a 
track length of 200 y. The average multiple-scattering angle is inversely 
proportional to the electrons momentum. For the case of sulfur, an average 
of 10 knock-on electrons/cm was found. This is compatible with the 
theoretical probability per unit pathlength of the production of a knock-on 
electron with energy e : 


o>(E,e ) 


2irNZ 2 q 4 

2 2 . ,2 
m c B (e ) 
e - 


S 2 (e_) 
^ e - ^max 


(1) 
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where N is the number of electrons per unit volume , E(0) is the incident 
energy of the primary ion, Z is the atomic number of the primary particle, and 
q is the electron charge. Also the emission angle is given by 


COS0 * 


8 T(e ) + 


2m c 
e 


( 2 ) 


and 


(e 


) 

- max 


M 


_ 2„2 2 
2m c g Y 
e 


(3) 


for m /M << 1; M is the mass of the heavy ion. Since y = 200, 

then fe ) ** 4.1 x 10 4 MeV, 6 “ 1 and so 

- max 


u(E,e ) 


2TrNZ 2 q 4 


m c 
e 



(4) 


Having found a knock-on electron, examine each origin grain for an associated 
track with e >1 MeV. This energy is a practical lower limit because the 
energy uncertainty approaches 50 percent. For a bonafide electron pair the 
origin grain of the electron (positron) must be within 1 y of the primary 
track and the origin grain of each associated track must be within 2 y of each 
other. With this criterion, it is estimated that 1-2 percent of the candidate 
pairs will be false. This background depends on the chance coincidence of a 
delta ray of energy e_ < 1 MeV and a fast knock-on electron of several MeV. 

A correlation plot (Figure 7) of the emission angle vs. the energy has 
been generated [27] for a subset (101 pairs) of the 200 GeV/amu sulfur 
electron pair data found with the latest scan method described above. Not 
unexpectedly, the maximum emission angle is inversely proportional to the 
energy but there is a large distribution of angles at each energy. The total 
energy transfer in the interaction (Figure 2), the energy division between the 
electron and positron (Figure 8), the net transverse momentum of the pairs 
(Figure 9), and the angular distribution for the electrons /positrons (Figure 
10) have all been derived from this same 200 GeV/amu sulfur electron pair data 
set. 


VI. DISCUSSION AND FUTURE WORK 

To resolve the question of whether the latest observed high yield of the 
low-energy pairs is contaminated by some background effect in the emulsion, 
improved calculations of the cross section for pair production by relativistic 
heavy nuclei [25] will be incorporated into a Monte Carlo simulation of the 
pair production process and compared to the data. The Monte Carlo simulation 
will be based on the Weizsacker-Williams virtual quanta method. The 
theoretical curves in Figures 2 and 10 are a result of the WW method. Further 
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measurements on oxygen ion tracks at 200 GeV/amu and 60 GeV/amu will be 
performed for this comparison. 


The differential cross sections that will be randomly sampled in the 
Monte Carlo simulation are 


da 


ww 

de 


N(e)*a (e), 

D 


(5) 


where N(e) is the Weizsacker-Williams frequency distribution for the virtual 
quanta, the energy transfer e = virtual photon energy, and a B (e) is the Racah 
expression [23] for the pair production by photons. 


d °NS _ Bethe-Heitler formula 
de + " 3D-1000 in [22]. 


This cross section, differential in positron energy, is based on the first 
Born approximation with a non-screened nucleus. 



f f 

J J de + dfi + 

e e , 


N(e )de + de 


(7) 


This cross section, differential in positron energy and solid angle, is the 
Sauter-Gluckstern-Hall formula (3D-2000 in ref. [22]) and is weighted with the 
Weizsacker-Williams distribution. Once again it applies for an unscreened 
point nucleus. 


A theory that has no free parameters and yet is applicable to the low- 
energy transfer pairs will be pursued. 

In the near future, an examination of some of the Japanese-American 
Collaborative Emulsion Experiment (JACEE) cosmic ray data will be done in 
order to compare the direct electron pair yield with the corresponding values 
calculated from the calorimeter energy estimates. Additional CERN experiments 
with a Pb beam, to be available in 1991, are planned since the increased 
electron pair yield may resolve any remaining observational problems as well 
as validate the calculations. 


This study has not yet elevated the measurement of the direct Coulomb 
pair yield to an energy measuring technique. It has, however, made 
considerable progress in both theory and observations that may allow in the 
near future the resolution of long standing discrepancies between various 
experimental results, and between calculations and observations of this 
interesting aspect of quantum electrodynamics. 
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TABLE 1. SUMMARY OF THE MEASURED DIRECT ELECTRON PAIR YIELD. (THE Fe 56 
DATA WERE TAKEN FROM THE LBL EXPOSURE WHILE THE O 1 6 AND S 32 DATA WERE 
ACQUIRED AT CERN. ) 


Projectile 

Primary Energy 
(GeV/nuc ) 

Track 

(cm) 

Pairs 

Jets 

Yield 

(pairs/cm) 

Fe 56 

1.8 

96.4 

27 

— 

0. 280±0. 050 

O 1 6 

60.0 

978.0 

46 

90 

0.047±0.007 

O 1 6 

200.0 

2959.0 

270 

218 

0.091 ±0.006 

s 32 

200.0 

979.6 

215 

120 

0.219±0.015 (a) 

s 32 

200.0 

129.4 

40 

16 

0.309±0.049 (b) 


(a) = 20 cm/day 

(b) = 5 cm/day 


9 





1000 


c 



NOisinwa ujoe/suivd N 


C 
03 
03 
42 

2 
<d 

>1 43 

P 03 


P 

e *— • 

<tS <Tk 

5 “ 


<d 

*§ 5 


O 

c 

£ 

03 

CM 

0 

■rH 

rH 

CO 


•P 

P 

03 0) 

id 

-P 

4-> 

04 

•P P 

o 

id 

O 


>, <d 



♦P 

CO 


O' 


T3 


0) CO 

p 

P 

• o 


>1 0) T3 
42 rH -P 
O rP 

p o 

C CO 


c 
o 

■p 

m >1 Q) 

2 g 


3 

£ 

03 

* 

o 

<d 

p 

4J 


03 

P 

a 

s 

s 

03 

C 

rP 

a; 


Li 03 

<d Li 

03 (U 

-h > un 

O 0 ) o 

3 CO 

C 6 
*0 U 
c c 

>p id co 


T> ^ C 

o> c *p 
U 0 
P P CO 
'O -H 3 
0 03 

Ll % rP 
04 P O 
P D 
(0 44 C 

Ll rH 
♦H P 03 
rd CO Cu 


C7» CO 03 

c <u ,c 

•P 4 -> 

c »p 

<D id o 

0) U -P 

Li *P 

0 4 - 3*0 

co a) a> 

P 43 
0 ) o u 

.5 a) -p 

5 5 i 

5 4,- 

O CO CN • 

42 o) cn • 

42 * VO 

Li * r— (N 
O <N 

mp a> •— ■ 

■g 

*- 43 

— . . i) O 

r— i 0) r- 

(d O E 

MW » 
id 

c 

c <d 


> 

x 


G 

O 

p 

V 

u 

0) 


o 

03 

u 

•P 

MH 

O 


e id 

0) 

tx» mh 
>i o 
X 

0 T3 

rH 

1 a> 

« -P 

§ * 

■P -P 

o c 

Ll 0> 

Qj rH 

<d 

* > 

CO 


id 

tH -P -rH 

I 43 *H 
Id ^ *-1 
•d 3 *H 
• Ll ^ 
■P I 
0 ) P 
PU QJ 
I M 
G O 


a) 


0) 

D 

3 

0 

Ll 

3 

£ 


C 

O 

•P 

O' 

a> 

Ll 

Ll 

03 


03 


C 

o 

Li 
-P 
^ O 

>i 0) 


a) a) 

fi w S 

CO 

03 03 
03 
N 
*P 


id mh 
CO CO 
PNC 
O *p id 
* 03 P 

03 

55 


03 MH 

_ P o 
•P *o 

XJ c 

03 CO O 
P *H *H 


>i 

Li 

0) 

C 

03 

03 

«P 

Id 

■P 


•H 
4 -> 

O 
03 Id 


03 

p _ _ 

p *o e 

O M 
•P P O 
Cm 04 C 


id 

o 

•P 
-P 
03 T3 
P 03 
O 43 
03 CO 


0) id T3 
> C 03 
6 


P . 

84 

O 


5 1.3 I 


10 




z a> 
< < 


11 


Figure 2. The energy transfer distribution for the direct pairs derived from a 
subset of the 200 GeV/amu Sulfur data [27]. The solid curve is a result of the 
Weizsacker-Williams-Kelner-Kotov calculation described in [26] . 
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Figure 5, A microphotograph mosaic of the 200 GeV/amu sulfur ion beam in emulsion 
including a direct electron pair originating on the left side of the photograph and 
traveling to the right. 
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DISTRIBUTION OF DIRECT ELECTRON PAIRS VS EMULSION DEPTH 
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Figure 6. Distribution of the direct electron pairs as a 
function of the distance from the edge of the emulsion plate 
for both the 200 GeV/amu oxygen and sulfur CERN exposures. 
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Figure 7 . A scatter plot of the individual emission angles of both the electron and 
positron of the direct pair versus their respective energies for a subset (101 
pairs) of the 200 GeV/amu sulfur data [27]. 
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Figure 9, The frequency distribution of the net transverse momentum for the direct 
electron pairs deduced from the subset of sulfur data referenced in Figure 7 • 
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